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Abstract 

In this talk, I present what I believe we have learned from the 
recent RHIC heavy ion experiments. The goal of these experiments is 
to make and study matter at very high energy densities, greater than 
an order of magnitude larger than that of nuclear matter. Have we 
made such matter? What have we learned about the properties of this 
matter? What do we hope and expect to learn in the future? 



1 Introduction 

The goal of the heavy ion program at the RHIC at Brookhaven National 
Laboratory is to make and study new forms of matter at energy densities in 
excess of ten times that of nuclear matter. I will describe the status of this 
program from a theorist's perspective. In this talk, I would like to address 
four simple questions: 

• What are we trying to understand? 



We are trying to produce new forms of matter and understand their 
properties. 

*Invited talk presented at the International Conference on the Physics and Astrophysics 
of the Quark Gluon Plasma, Nov. 26-30 2001 
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• What have we already learned? 

I will argue that we have produced matter at energy densities at least 
10 and perhaps 100 times higher than that of nuclear matter. This 
matter is strongly interacting with itself. 

• What do we expect to learn? 

There are some measurements which will be carried out in the near 
term which will answer specific questions about the properties of matter 
produced in RHIC collisions. 

• What do we hope to learn? 

In the longer term, there are the measurements which are harder to 
understand or more controversial or ambiguous in their interpretation, 
and will require new analysis. 

2 What Are We Trying to Understand? 

There are two central issues of the RHIC experimental program. 

• What is the behavior of matter at asymptotically large energy 
density? 

Matter at very high energy density in thermal equilibrium is believed 
to form a Quark Gluon Plasma. This gas of almost free quarks and 
gluons is thought to be the proper description of matter when energy 
densities are larger than about 1 GeV/Fm 3 . This is about the energy 
density inside a proton or neutron and is about an order of magnitude 
larger than that of nuclear matter. Matter of this energy density occurs 
naturally in the cores of neutron stars and was present during the big 
bang. 

• What type of matter is important for high energy hadrons? 
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This matter is believed to be a very dense system of quarks and glu- 
ons, and is called a Color Glass Condensate. It is an incoherent su- 
perposition of quantum mechanical states of Bose condensates. It is 
parameterized by a surface energy density, and is presumably the cor- 
rect description of matter when the energy density per unit area is 
greater than about 1 GeV/Fm 2 The Color Glass Condensate may be 
responsible for universal behavior of all hadronic interactions at high 
energy. 



2.1 The Quark Gluon Plasma 

In Fig. [I], a cartoon is shown which illustrates what we expect happens when 
nucleons, or for that matter any hadrons, are compressed to densities higher 
than is typical inside a nucleon. Eventually the constituents of the nucleon 
travel more or less freely around the system as a whole. This system is the 
Quark Gluon Plasma. 




Compressed Nucleons form Quark Gluon Plasma 



Figure 1: Nuclear matter when compressed becomes a gas of quarks and 
gluons: the Quark Gluon Plasma 
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The formation of a Quark Gluon Plasma should begin to occur at about 
the energy density of matter inside a proton. Since a proton has a size about 
R ~ 1 Fm, and has a rest mass energy of M ~ 1 GeV, corresponding to an 
energy density of about 1 GeV/Fm 3 

These observations led a number of people in the late 1970's to suggest 
that there was a phase transition between ordinary matter and a Quark 
Gluon Plasma as is shown in Fig. |2|.[|l[] These conjectures were later firmed 
up by lattice Monte-Carlo simulations. 0] 



The Evolving QCD Phase Transition 

t ~ 1980 

Quark Gluon Critical Temperature 150 - 200 MeV ( H b = ) 
Plasma Critical Density 1/2-2 Baryons/Fm 3 ( T = ) 
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Figure 2: Beliefs about the QCD phase diagram as a function of time. Solid 
line indicates phase boundary, and dashed a rapid crossover. 



The early lattice Monte-Carlo simulations were unable to properly in- 
clude the effects of quark masses. About the mid 1980's arguments ap- 
peared which strongly suggested that for realistic values of quark masses, a 
real phase transition does not occur. Instead there is a sharp change in the 
properties of the system at some temperature and baryon chemical poten- 
tial. Although it does not have strong implications for heavy ion collisions, 
this would affect strongly cosmological scenarios and has potential effect in 
neutron star cores. 
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By 2000, the story had become more complicated. It is now believed 
that there is a line of first order phase transitions in the baryon number 
chemical potential and temperature plane. There may be some region where 
the transition is second order or crossover. At high baryon number density 
and small temperature, there may be a number of color superconducting 
phases. 

The bottom line on this historical aside is that physics is an experimental 
discipline. I think it is unlikely we will have a first principles understanding 
of the phase transition in QCD on the basis of pure thought. Much work will 
have to be done which combines experimental results and difficult numerical, 
presumably lattice, Monte-Carlo simulation, before we will have a compelling 
picture of high density matter. 

At zero baryon number density, a lot has been learned about the proper- 
ties of high temperature matter. For technical reasons, lattice Monte Carlo 
methods are very difficult to implement at finite baryon number density. At 
finite temperature, one has made remarkable progress. || 

Massive Nucleoli 




Figure 3: (a) The energy density vs temperature, (b) The chiral condensate 
vs T. Both plots are at zero baryon number density. 

In Fig. |^a, a plot of the energy density, e scaled by T 4 is shown vs 
temperature, T. This variable e/T 4 is roughly proportional to the number of 
particle degrees of freedom important at the energy scale corresponding to 
T. A rapid cross over occurs at a temperature of T ~ 160 — 190 MeV. A first 
order phase transition cannot be ruled out at this time because the effects 
of finite quark mass are not under sufficient control. At low temperatures, 
the number of degrees of freedom are small and consistent with a gas of 
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pions. At high temperatures, the degrees of freedom are roughly that of a 
massless gas of quarks and gluons. A sharp transition occurs between the 
meson system and that of the quark and gluons where e/T 4 changes by an 
order of magnitude. By a temperature of order T ~ 500 MeV, the number 
of important particle degrees of freedom stays approximately constant. 

The origin of the proton and neutron mass is one of the mysteries of 
QCD. The masses of the up and down quarks inside a proton are only a few 
percent of the nucleon mass. Do not be fooled by those who say that the 
LHC Higgs boson search is designed to explain the origin of mass. The LHC 
probes electroweak physics and this gives masses to the up and down quark 
masses, which can explain only a few percent of the nucleon mass. Even in 
a world where the up and down quark masses had no mass and electroweak 
physics was entirely ignored, the nucleon would still have roughly its present 
mass, and this arises from QCD. It is believed that this mass is obtained by 
breaking of a chiral symmetry of the strong interactions. 




Figure 4: (a) The potential between two quarks as a function of distance. 
(b)The sound velocity as a function of temperature. 

This symmetry is broken in our low temperature world but is restored 
in the high temperature world. A measure of this symmetry breaking is 
< >, the typical value of a condensate. This condensate is composed 
of quark-antiquark pairs. At low temperatures, this is non-zero, and rapidly 
goes to zero at the phase transition temperature, as is seen in Fig. |d. The 
nucleon mass is proportional to this condensate, and so it goes to zero in the 
Quark Gluon Plasma. 

Showing that the confining force between two quarks is linear was one 
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of the first triumphs of lattice QCD. It can be measured as a function of 
temperature, and indeed in the Quark Gluon Plasma it goes to a constant 
at long distances. This is shown in Fig. [|a. 

How is this transition from a gas of pions into a Quark Gluon Plasma 
manifest in terms of bulk properties of the system? The sound velocity is 
v% = dP/de where P is the pressure. At the crossover, the pressure is roughly 
constant while the energy density changes by an order of magnitude. The 
sound velocity should drop to near zero at T c , as shown in Fig. f|b. 

2.2 The Color Glass Condensate 




To understand this new form of matter, it is convenient to imagine a 
hadron in a reference frame where it has very large longitudinal momentum. 
We will be interested in the constituents of the hadron wavefunction which 
have small longitudinal momentum in this frame of reference. These low mo- 
mentum constituents are produced by the high momentum ones. Because the 
high momentum constituents appear to have time scales which are Lorentz 
time dilated compared to their natural scales, and since they induce the low 
momentum fields associated with the low momentum particles, the low mo- 
mentum fields evolve very slowly compared to their natural time scale. Hence 
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the term, Color Glass, since the fields are composed of color gluons, and glass 
because the time scale for evolution of these low momentum fields is much 
longer than their natural time scale. These fields live on a two dimensional 
sheet because of the Lorentz contraction of the high energy hadron. We shall 
argue in the following paragraphs that the phase space-density of these fields 
becomes large and forms a condensate. ||] 

The fields on the two dimensional sheet turn out to be similar to the 
Lienard-Wiechart potentials of electrodynamics. They correspond to plane 
waves as in the Weizsacker- Williams approximation of electrodynamics, ex- 
cept that they have color. They have their color electric field perpendicular 
to their color magnetic field and both perpendicular to their direction of 
motion, E a _L B a _L z. They have a random color. This is shown in Fig. |5| 

The gluon structure function xG(x, Q 2 ) is experimentally measured to 
increase at small x. In the reference frame where the hadron is very fast, 
x is the ratio of a constituent energy to the projectile energy. The gluon 
distribution is shown in Fig. |6]a. Note the rapid increase in xG(x, Q 2 ) as a 
function of x for small x. This is the origin of the "small x problem" . This 
means that the piece of the hadron wavefunction relevant for small x processes 
has an increasing density of gluons. In Fig. |6|b, we look at a hadron headed 
along the beam direction. As x decreases, the density of gluons increases. 




Figure 6: (a) The gluon structure function as a function of x for various Q 2 . 
(b)The increase in density of gluons as x decreases. 
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The phase space density of gluons is 



1 dN 



(1) 



P 



ttR 2 dyd 2 pT 



where R is the hadron size, px is the transverse momentum of a constituent, 
and y ~ ln{l/x). The high density of gluons is generated dynamically and is 
caused by an instability, which is proportional to the density. The instability 
is stabilized when the density of partons becomes large enough so that inter- 
actions of order giqcdT) 2 become of the order of the linear instability Here 
rj — J d 2 pxP- This requires that 



The factor of Q 2 arises because we consider densities per unit area, and Q 2 at 
carries this dimension. This Q sat is called the saturation momentum. The 
factor of olqcd is the strong coupling strength of QCD. When Q sat » Aq C d, 
we expect that olqcd « 1, so that the system becomes a high density Bose 
Condensate. 

The name Color Glass Condensate arises therefore because 

• Color 

The gluons are colored. 

• Glass The natural time scale for the evolution of the gluon field is 
Lorentz time dilated. This is like a glass which is a liquid on long times 
scales but a solid on short ones. 

• Condensate 

The phase space density is as large as it can be. 

2.3 Space Time Evolution of Heavy Ion Collisions 

A collision of two sheets of Colored Glass is shown in Fig. 

This is the picture of nucleus-nucleus collision which arises from the Color 
Glass Condensate. 

The time evolution of the matter produced in these collisions is divided 
into several stages: 



V ~ Qlat/ a QCD 



(2) 
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Figure 7: High energy nucleus-nucleus collisions. 



• Initial Conditions 

For t < 0, the two sheets approach one another. The Color Glass is 
frozen in each nucleus. 

• Melting the Color Glass 

During the time < t < tf orm , the Color Glass melts into quarks and 
gluons. It is estimated that tf orm ~ 1/Q sa t ~ .1 — .3 Fm/c at RHIC 
energy. The energy density of the matter at formation is somewhere 
around €f orm ~ Q 4 sat /a s ~ 20 - 100 GeV/Fm 3 . 

• Thermalization 

During the time tf orm < t < t therm , the matter expands and ther- 
malizes. Typical thermalization time is estimated to be ttherm ~ -5 — 
1 Fm/c. 

• Hydrodynamic Expansion 

The system expands as a thermal system until a time of decoupling 
which is typically about t decoupling ~ 10 Fmj c at RHIC energy. Here the 
matter presumably starts as a Quark Gluon Plasma, evolves through 
a mixed phase of hadrons and Quark Gluon Plasma and eventually 
becomes a gas of pions. In this stage, most of the physics interesting 
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for studies of the phase transition or cross over between Quark Gluon 
Plasma and ordinary hadronic matter takes place. 

As the Color Glass melts, it produces particles as is shown in Fig. |[ 
The fastest particles have their natural time scale time dilated the most, so 




large p small p large p 



Figure 8: Particle production in nucleus- nucleus collision in the center of 
mass frame. 

in the center of mass frame, the fastest particles are produced last. These 
particles have traveled the longest distance from the collision point, since 
their formation time is dilated the most. Therefore the matter is formed 
with a correlation between momentum and position. This is like Hubble flow 
in cosmology where if you look at stars, the stars which are farthest away 
move away the fastest. For heavy ions as in cosmology, this description is 
frame independent. Unlike in cosmology, the Hubble expansion for heavy ion 
collisions is 1 dimensional. 0] 

The density of particles falls as N/V ~ 1/t. If the particles expand 
without interaction, then the energy per particle is constant. If the particles 
thermalize, then E/N ~ T, and since N/V ~ T 3 for a massless gas, the 
temperature falls as T ~ t -1 / 3 . For a gas which is not quite massless, the 
temperature falls somewhere in the range T Q > T > T (t /ty/ 3 , that is the 
temperature is bracketed by the value corresponding to no interaction and 
to that of a massless relativistic gas. This 1 dimensional expansion continues 
until the system begins to feel the effects of finite size in the transverse 
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direction, and then rapidly cools through three dimensional expansion. Very 
close to when three dimensional expansion begins, the system decouples and 
particle free stream without further interaction to detectors. 



3 What Have We Learned? 
3.1 The Energy Density is Big 

The particle multiplicity as a function of energy has been measured at RHIC 
and is shown in Fig. [)|[|9[]-|12|] On the same plot is shown lower energy data 
from the AGS and CERN and data from the pp collider as well. One sees 
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Figure 9: The multiplicity as a function of energy for heavy ions and for pp. 



that the energy dependence in AA collisions is different than that for pp, and 
is more or less consistent with a ln(s) behavior. 

Combining the multiplicity data together with the measurements of trans- 
verse energy or of typical particle transverse momenta, one can determine the 
energy density of the matter when it decouples. One can then extrapolate 
backwards in time. We extrapolate backwards using 1 dimensional expan- 
sion, since decoupling occurs when the matter first begins to expands three 
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Figure 10: Bounds on the energy density as a function of time in heavy ion 
collisions. 



dimensionally. We can extrapolate backwards until the matter has melted 
from a Color Glass. We shall take a conservative overestimate of this time 
to be of order t me i t ~ .3 Fm/c The extrapolation backwards is bounded by 
€f(tf/t) < e(t) < e/(t//t) 4 / 3 . The lower bound is that assuming that the 
particles do not thermalize and their typical energy is frozen. The upper 
bound assumes that the system thermalizes as an ideal massless gas. We 
argued above that the true result is somewhere in between. When the time 
is as small as the melting time, then the energy density begins to decrease 
as time is further decreased. 



This bound on the energy density is shown in Fig. [10]. On the left axis is 
the energy density and on the bottom axis is time. The system begins as a 
Color Glass Condensate, then melts to Quark Gluon Matter which eventually 
thermalizes to a Quark Gluon Plasma. At a time of a few Fm/c, the plasma 
becomes a mixture of quarks, gluons and hadrons which expand together. At 
a time of about 10 Fm/c, the system falls apart and decouples. At a time of 
t ~ 1 Fm/c, the estimate we make is identical to the Bjorken energy density 
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Figure 11: The multiplicity per participant as a function of the number of 
participants. The solid line running through the center of both of the sets 
of data is the Kharzeev-Nardi result. The EKRT result is marked as such in 
the left hand plot and is called the saturation model in the right hand plot. 
HIJING is marked on both plots. 



estimate, and this provides a lower bound on the energy density achieved in 
the collision. (All estimates agree that by a time of order 1 Fm/c, matter has 
been formed.) The upper bound corresponds to assuming that the system 
expands as a massless thermal gas from a melting time of .3 Fm/c. (If the 
time was reduced, the upper bound would be increased yet further.) The 
bounds on the energy density are therefore 

2-3 GeV/Fm 3 < e < 20 - 100 GeV/Fm 3 (3) 

where we included a greater range of uncertainty in the upper limit because 
of the uncertainty associated with the formation time. The energy density of 
nuclear matter is about 0.15 GeV/Fm 3 , and even the lowest energy densities 
in these collisions is in excess of this. At late times, the energy density is 
about that of the cores of neutron stars, e ~ 1 GeV/Fm 3 . 

At such extremely high energy densities, it is silly to try do 
describe the matter in terms of anything but its quark and gluon 
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degrees of freedom. 



3.2 The Gross Features of Multiplicity Distributions 
Are Consistent with Colored Glass 
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Figure 12: Color glass condensate fits to the rapidity density measured in 
the PHOBOS and Brahms experiments 



As argued by Kharzeev and Nardi,|T3[ the density of produced particles 
per nucleon which participates in the collision, N part , should be proportional 
to l/a s (Qsat), an d Q 2 sat ~ N part . This dependence follows from the l/a s 
which characterizes the density of the Color Glass Condensate. In Fig. |ll], we 
show the data from PHENIX and PHOBOS [|14]]. The Kharzeev-Nardi form 



fits the data well. Other attempts such as HIJING[]T5[, or the so called satu- 



ration model of Eskola-Kajantie-Ruuskanen-Tuominen|l6[ are also shown in 
the figure. 

Kharzeev and Levin have recently argued that the rapidity distribu- 
tions as a function of centrality can be computed from the Color Glass 
description. [O] This is shown in Fig. O. 
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3.3 Matter Has Been Produced which Interacts Strongly 
with Itself 



In off zero impact parameter heavy ion collisions, the matter which overlaps 
has an asymmetry in density relative to the reaction plane. This is shown in 
the left hand side of Fig. [TB|. Here the reaction plane is along the x axis. In 
the region of overlap there is an x — y asymmetry in the density of matter 
which overlaps. This means that there will be an asymmetry in the spatial 
gradients which will eventually transmute itself into an asymmetry in the 
momentum space distribution of particles, as shown in the right hand side 



of Fig. 13. This asymmetry is called elliptic flow and is quantified by the 




In region of overlap 
asymmetry of spatial 
x distribution -> 
asymmetry of momentum 
distribution 



V, = < cos(2 <|> )> tan( § ) = Py/p x 



Figure 13: The asymmetry in the distribution of matter in an off center 
collision is converted to an asymmetry of the momentum space distribution. 



variable t> 2 . in all attempts to theoretically describe this effect, one needs 
very strong interactions among the quarks and gluons at very early times in 
the collision. |19| . In Fig. 14], two different theoretical descriptions are fit to 



the data by STAR and PHOBOS|3-0. On the left hand side, a hydro- 
dynamical model is used.|22|] It is roughly of the correct order of magnitude 
and has roughly the correct shape to fit the data. This was not the case 
at lower energy. On the right hand side are preliminary fits assuming Color 
Glass. [23] Again it has roughly the correct shape and magnitude to describe 
the data. In the Color Glass, the interactions are very strong essentially from 
t = 0, but unlike the hydrodynamic models it is field pressure rather than 
particle pressure which converts the spatial anisotropy into a momentum 
space-anisotropy. 
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Probably the correct story for describing flow is complicated and will in- 
volve both the Quark Gluon Plasma and the Color Glass Condensate. Either 
description requires that matter be produced in the collisions and that it in- 
teracts strongly with itself. In the limit of single scatterings for the partons 
in the two nuclei, no flow is generated. 
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Figure 14: a) A hydrodynamic fit to v2. b) The Colored Glass fit. 



4 What Do We Expect to Learn? 
4.1 Does the Matter Equilibrate? 

One of the most interesting results from the RHIC experiments is the so 
called "jet quenching" . [plf l- 125| . In Fig. [l5la, the single particle hadron spec- 
trum is scaled by the same result in pp collisions and scaled by the number 
of collisions. The number of collisions is the number of nucleon-nucleon in- 
teractions, which for central collisions should be almost all of the nucleons. 
One is assuming hard scattering in computing this number, so that a single 
nucleon can hard scatter a number of times as it penetrates the other nucleus. 
The striking feature of this plot is that the ratio does not approach one at 
large px- This would be the value if these particles arose from hard scattering 
which produced jets of quarks and gluons, and the jets subsequently decayed. 

The popular explanation for this is shown in Fig. |T5|b. Here a pair of 
jets is produced in a gluon-gluon collision. The jets are immersed in a Quark 
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Figure 15: a) The px distribution of particles scaled by the data from pp col- 
lisions times the number of hard collisions inside the nuclei, b) A pair of jets 
is produced in a hard collision and they propagate through the surrounding 
matter. 

Gluon Plasma, and rescatter as they poke through the plasma. This shifts 
the transverse momentum spectrum down, and the ratio to pp collisions, 
where there is no significant surrounding media, is reduced. 

The data, however suggestive, need to be improved before strong conclu- 
sions are drawn. For example, there are large systematic uncertainties in the 
pp data which was measured in different detectors and extrapolated to RHIC 
energy. This will be resolved by measuring pp collisions at RHIC. There is in 
addition some significant uncertainty in the AA data which becomes smaller 
in the ratio to pp data when the data is measured in the same detector. 
There are nuclear modifications of the gluon distribution function, an effect 
which can be determined by measurements on pA at RHIC. The maximum 
transverse momentum is limited by the event sample size, and the size will 
be greatly improved with this years run due to the higher luminosity and 
longer run time. 

One of the reasons why jet quenching is so important for the RHIC 
program is that it gives a good measure of the degree of thermalization in 
the collisions. If jets are strongly quenched by transverse momenta of 4 GeV, 
then because cross sections go like 1/E 2 for quarks and gluons, this would be 
strong evidence for thermalization at the lower energies typical of the emitted 
particles. 

One can look for evidence of thermalization directly from the measured 
Pt distributions. Here one can do a hydrodynamic computation and in so 
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far as it agrees with the results, one is encouraged to believe that there is 
thermalization. On the other hand, these distributions may have their origin 
in the initial conditions for the collision, the Colored Glass. In reality, one 
will have to understand the tradeoff between both effects. 

Let us begin with the Color Glass description. In Fig. |l6|a, the mT 
distributions (m^ = p\ + m 2 ) of identified particle measured for minimum 
bias events in PHENIX are shown. |25[] In Fig [T6|b, we replot these curves 
by rescaling by a constant multiplicative factor for each of the particles. 
In thermal models, this constant factor would be associated with chemical 
potentials for producing different particle species. These plots fall on top of 
one another. The distributions therefore scale in tut- 
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Figure 16: a The itlt distributions of identified particles, b) The distributions 
rescaled by a constant dependent upon particle species. 



If the distributions scale in tut, as they appear to do for itlt within 
the range measured by PHENIX, then the Colored Glass predicts as scaling 
relation for the distributions: 



dN 
dycPpr 



ttR* 



F(m T /Q S at) 



(4) 



This means that if we adjust the normalization of the distribution function 
and the scale of its dependence on ttit then the distributions should fall on 
top of one another. The scale adjustment should be proportional to l/a s 
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which should in turn involve Q sa t from the running of the coupling. In Fig. 
17 a, the data at different centralities are rescaled into one another according 
to the above relation. This seems to describe the data well. In Fig. [TTpo, the 
dependence of l/a s extracted from the above equation is compared to what 
is determined from the running coupling constant. Again, this works quite 
well. The PHENIX data is consistent with the scaling relations predicted by 
the Color Glass Condensate. 
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Figure 17: a) The ray distributions for various centralities have been rescaled 
into one another, b) The dependence of the coupling l/a s is compared to 
the empirically determined scaling factor for the rrij- distributions. 



Hydrodynamical models also successfully describe the data on m-r 
distributions. |27j In Fig. [H| the results of the simulation by Shuryak and 
Teaney are shown compared to the STAR and PHENIX data.p^|-|[25|| The 
shape of the curve is a prediction of the hydrodynamic model, and is param- 
eterized somewhat by the nature of the equation of state. Notice that the 
STAR data include protons near threshold, and here the tut scaling breaks 
down. The hydrodynamic fits get this dependence correctly, and this is one 
of the best tests of this description. The hydrodynamic models do less well 
on fits to the more peripheral collisions, but the Colored Glass model gets 
this more or less correctly. In general, a good place to test the hydrodynamic 
models predictions is with massive particles close to threshold, since here one 
deviates in a computable way from the rriT scaling curve, which is arguably 
determined from parameterizing the equation of state. 
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Figure 18: The hydrodynamical model fits to the mx spectra for the PHENIX 
and STAR data. 



If one can successfully argue that there is thermalization in the RHIC 
collisions, then the hydrodynamic computations become compelling. One 
should remember that hydrodynamics requires an equation of state plus ini- 
tial conditions, and these initial conditions are determined by Colored Glass. 
Presumably, a correct description will require the inclusion of both types of 
effects. |2§ 

5 What Do We Hope to Learn? 

5.1 Confinement and Chiral Symmetry Restoration 

We would like to know whether or not deconfinement has occurred in dense 
matter or whether chiral symmetry restoration has changed particle masses. 

This can be studied in principle by measuring the spectrum of dileptons 
emitted from the heavy ion collision. These particles probe the interior of the 
hot matter since electromagnetically interacting particles are not significantly 
attenuated by the hadronic matter. For electron-positron pairs, the mass 



21 




Figure 19: The CERES data on low mass electron-positron pairs. The ex- 
pected contribution from ordinary hadrons is shown by the solid line. The 
data points are for the measured electron-positron pairs. 



distribution has been measured in the CERES experiment at CERNf29||, and 
is shown in Fig. [19]. Shown in the plot is the distribution predicted from 
extrapolating from pA collisions. There should be a clear p and <p peak, which 



has disappeared. This has been interpreted as a resonance mass shift, [30 



enhanced rj' production, but is most probably collisional broadening of 
the resonances in the matter produced in the collisions. [32] Nevertheless, if 
one makes a plot such as this and the energy density is very high and there 
are no resonances at all, then this would be strong evidence that the matter 
is not hadronic, i. e. the hadrons have melted. 

The resolution in the CERES experiment is unpleasantly large, making 
it difficult to unambiguously interpret the result. Whether or not such an 
experiment could be successfully run at RHIC, much less whether the reso- 
lution could be improved, is the subject of much internal debate among the 



22 



RHIC experimentalists. 



5.2 Confinement and J/ty Suppression 



In Fig. |2(J, the NA(50) data for J/\& production is shown. [33] In the first fig- 
ure, the ratio of J/\I/ production cross section to that of Drell-Yan is shown as 
a function of E T , the transverse energy, for the lead-lead collisions at CERN. 




Figure 20: a) The ratio of produced J pairs to Drell-Yan pairs as a function 
of transverse energy Et at CERN energy, b) The measured compared to the 
theoretically expected J/\l/ suppression as a function of the Bjorken energy 
density for various targets and projectiles. 



There is a clear suppression at large which is greater than the hadronic 
absorption model calculations which are plotted with the data.[|4]] In the 
next figure, the theoretically expected J/\l/ suppression based on hadronic 
absorption models is compared to that measured as a function of the Bjorken 
energy density for various targets and projectiles. There appears to be a 
sharp increase in the amount of suppression for central lead-lead collisions. 

Is this suppression due to Debye screening of the confinement potential in 
a high density Quark Gluon Plasma? ||35||-|37|| Might it be due higher twists, 



enhanced rescattering, or changes in the gluon distribution function? [|38 



[3"9|] Might the J/ip suppression be changed into an enhancement at RHIC 
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energies and result from the recombination in the produced charm particles, 
many more of which are produced at RHIC energy? [40|j-[ 43[ 



These various descriptions can be tested by using the capability at RHIC 
to do pp and pA as well as A A. Issues related to multiple scattering, higher 
twist effects, and changes in the gluon distribution function can be explored. 
A direct measurement of open charm will be important if final state recom- 
bination of the produced open charm makes a significant amount of J/^'s. 

5.3 The Lifetime and Size of the Matter Produced 

The measurement of correlated pion pairs, the so called HBT pion inter- 
ferometry, can measure properties of the space-time volume from which the 



hadronic matter emerges in heavy ion collisions. j44J The quantities Ri on g, Rside 
and Rout shown in Fig. |21] measure the transverse size of the matter at de- 
coupling and the decoupling time. 




K out 

Figure 21: The various radii used for HBT pion interferometry. 



In Fig. [221, the data from STAR and PHENIX is shown. |45||- [}4q| There 
is only a weak dependence on energy, and the results seem to be more or 
less what was observed at CERN energies. This is a surprise, since a longer 
time for decoupling is expected at RHIC. Perhaps the most surprising result 
is that Rout/ Rside is close to 1, where most theoretical expectations give a 
value of about R ou t/ Rside ~ 2.j47|-f48| Perhaps this is due to greater than 
expected opacity of the emitting matter? At this time, there is no consistent 
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Figure 22: a) The various HBT radii measured in heavy ion experiments 
including the new data from STAR, b) The correlation functions which de- 
termine the radii as a function of the pair momenta measured in PHENIX. 



theoretical description of the HBT data at RHIC. Is there something missing 
in our space-time picture? 

5.4 The Flavor Composition of the Quark Gluon Plasma 

The first signal proposed for the existence of a Quark Gluon Plasma in heavy 
ion collisions was enhanced strangeness production. |49| This has lead to a 



comprehensive program in heavy ion collisions to measure the ratios of abun- 
dances of various flavors of particles. [EO] . In Fig. p3]a, the ratios of flavor 



abundances is compared to a thermal model for the particle abundances. [51 



The fit is quite good. In Fig. |2"3|b, the temperature and baryon chem- 
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Figure 23: a) Various ratios of particle abundances and the RHIC data. The 
lines are the predictions of a thermal model, b) The temperature vs baryon 
chemical potential for a thermal model which is fit to data at various energies. 
The dashed line is a hypothetical phase boundary between a Quark Gluon 
Plasma and a hadronic gas. 



ical potential extracted from these fits is shown for experiments at various 
energies and with various types of nuclei. It seems to agree reasonably well 
with what might be expected for a phase boundary between hadronic matter 
and a Quark Gluon Plasma. 

This would appear to be a compelling case for the production of a Quark 
Gluon Plasma. The problem is that the fits done for heavy ions to particle 
abundances work even better in e + e~ collisions. One definitely expects no 
Quark Gluon Plasma in e + e~ collisions. There is a deep theoretical question 
to be understood here: How can thermal models work so well for non-thermal 
systems? Is there some simple saturation of phase space? The thermal model 
description can eventually be made compelling for heavy ion collisions once 
the degree of thermalization in these collisions is understood. 



6 Summary 

• What are we trying to understand? 

We want to make and study new forms of high energy density matter. 
At early times in a heavy ion collisions, this might be a Color Glass 
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Condensate, and at later times a Quark Gluon Plasma. 

• What have we already learned? 

Matter has been produced at energy densities so high it can only be sim- 
ply described in terms of quarks and gluons. This matter is strongly in- 
teracting with itself. The multiplicity distributions of matter produced 
in heavy ion collisions is consistent with it being a Colored Glass. Flow 
and pt distributions are consistent with both a Quark Gluon Plasma 
and a Colored Glass. 

• What do we expect to learn? 

By measuring the pr distributions of heavy hadrons, we can resolve the 
differences between various theoretical descriptions of the evolution of 
the produced matter. By measuring high pt hadrons in AA and pA 
collisions, we can determine whether jet quenching occurs, and to what 
degree the matter has thermalized. 

If thermalization has occurred, we can then perhaps measure the equa- 
tion of state. This is important because the equation of state reflects 
the bulk properties of the matter, and determines hydrodynamic and 
static properties of Quark Gluon Matter. 

• What do we hope to learn? 

Through measurements of lepton pairs, we can determine whether 
or not resonances have shifted in mass, have broadened or melted. 
Through measurements of pion pairs, we can find the spatial size and 
lifetime of the matter produced in heavy ion collisions, and test space- 
time pictures of ultrarelativistic nuclear collisions. Flavor abundances 
can tell us about the macroscopic properties of hadronic matter. The 
J/ty abundance may tell us whether the matter produced in these col- 
lisions is confining or not. 

We must turn to the problem of determining the properties of 
matter made in heavy ion collisions. We have made matter at such a 
high energy density that it can only be simply described in terms of quarks 
and gluons. The future will be to understand the properties of this matter. 

At present we are developing the theoretical and experimental tools to 
analyze ultra-relativistic heavy ion collisions. It is a learning time for both 
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theorists and experimentalists. Ultimately, we will need experiments with pp, 
pA and AA to resolve the initial state effects due to a Color Glass Condensate 
from effects related to the evolution of the matter as a Quark Gluon Plasma. 
A more precise determination of the properties of Color Glass Condensate 
would be possible with a high energy eA collider or in the LHC heavy ion 
experiments. 
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